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摘要 : 嗜 尸 性 蝇 类 在 命案 死亡 时 间 和 现场 推断 方面 有 着 十 分 重要 的 应 用 , 而 DNA. 条 形 编 码 技术 能 摆脱 对 忠 卵 和 幼 
忠 的 饲养 以 及 后 续 物 种 鉴定 方面 专业 知识 的 依赖 , 有 助 于 实现 现场 采集 蝇 类 样本 的 快速 鉴定 。 本 人 研究 采集 了 北京 地 
区 7 个 嗜 尸 性 蝇 类 优势 种 共 TT 个 个 体 的 样本 , 测定 了 所 有 个 体 线粒体 DNA. 上 细胞 色素 C 氧化 酶 亚 基 工 (COI) 基因 
1 120 bp 的 序列 。 基 于 序列 的 系统 发 生 分 析 显 示 , 同一 物种 不 同 个 体 的 序列 均 以 高 达 99% 的 文 持 值 聚 集 在 一 起 。 序 
列 间 的 分 长 统计 表明 这 些 蝇 类 在 物种 内 的 个 体 分 歧 不 超过 196 , 而 不 同 物种 间 的 净 分 歧 均 超过 7.74% , 最 高 可 达 
14.8596 。 滑 动 窗口 分 析 表 明 , 在 整个 序列 区 段 种 间 差 异 位 点 存在 较 平均 的 分 布 。 通 过 测定 COI 基因 的 序列 ,建立 
了 北京 地 区 7 个 嗜 尸 性 蝇 类 优势 种 的 DNA 条 形 码 , 据 此 实现 了 对 这 些 物 种 准确 快速、 简单 的 区 分 和 鉴定 , 同时 也 为 
后 续 应 用 于 物种 鉴定 的 种 属 特 异性 位 点 之 筛选 提供 了 基础 数据 。 
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Abstract: Sarcosaphagous flies have very important applications in death time estimation and crime scene 
reconstruction for homicide cases. DNA barcoding technology can extricate the investigators from the 
difficulty of culturing eggs and larvae of these flies as well as the reliance on special knowledge to 
scientifically identify the following imagoes which is of great help for rapid identification of fly samples 
collected at the crime scene. We collected 77 samples, which belong to 7 dominant species of 
sarcosaphagous flies in Beijing, from different districts, and sequenced a 1 120 bp region of the 
mitochondrial cytochrome oxidase subunit I ( COI) gene for each individual. Phylogenetic analysis showed 
that individual sequences of each species highly clustered together with 99% bootstrap support. Sequence 
divergences within species were less than 1% , while net divergences between species were all greater than 
7.7496 , the highest reaching 14. 8596. Sliding window analysis demonstrated that variable sites for all 


species were relatively equally distributed along the entire region. By sequencing part of the COI gene, we 
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established the DNA barcodes for the seven dominant species of sarcosaphagous flies in. Beijing, which 


realized the accurate, rapid, and easy identification of these species and provided basic material for the 


future selection of variable sites that can be used for direct identification of sarcosaphagous flies. 


Key words: Forensic entomology; sarcosaphagous flies; species identification; mitochondrial COI gene; 


DNA barcode; Beijing 


1996 年 北京 市 公安 局 与 中 国 科 学 院 动物 研究 
所 合作 , 对 北京 地 区 哮 尸 性 晶 类 展开 了 系统 的 研 
究 , 在 几 年 内 连续 采集 了 北京 不 同 地 区 、 不 同 季 市 
( 春 夏 秋 ) 和 不 同 环境 下 命案 尸体 周围 出 现 的 嗜 己 
性 蝇 类 的 不 同 虫 态 , 通过 对 成 虫 的 固定 、 幼 虫 及 虫 
卵 的 饲养 并 进行 最 终 的 形态 学 鉴定 ， 总结 出 北京 地 
区 哮 尸 性 蝇 类 优势 种 涉及 3 科 6 属 共 7 RR CE SER 
等 ,1998) : 红头 丽 蝇 Calliphora vicina, 大头 金 量 
Chrysomyia megacephala 丝光 绿 量 Lucilia sericata HE 
2l MV. fk te Parasarcophaga crassipalpis . 3 £4] MV. EK 98 
Parasarcophaga portschinskyi Rohdendorf. f Fé 5!) Jk tE 
Boettcherisca. peregrina 和 家 晶 Musca domestica ~ 

传统 物种 鉴定 工作 主要 依靠 形态 学 上 的 特征 ， 
但 是 用 于 区 分 物种 的 形态 学 差别 很 多 时 候 是 如 此 之 
复杂 和 细微 , 以 至 于 经 常 仅 有 少数 专家 能 对 一 组 关 
系 紧密 的 物种 进行 准确 区 分 。 还 有 一 些 物种 更 多 以 
幼 态 (juvenile form) 形式 存在 ， 而 此 时 它们 并 不 表 
现 出 任何 物种 水 平 的 区 分 特征 , 只 有 等 培养 成 熟 以 
后 (如 果 可 以 的 话 ) 才 能 被 准确 鉴定 (Catts and Goff, 
1992; Stoeckle, 2003) 。 作 为 应 用 学 科 的 法 医 昆 虫 
学 , 现场 发 现 的 蝇 卵 或 蛆虫 往往 要 人 饲养 一 段 时 间 后 
才能 进行 形态 学 种 属 鉴定 , 这 对 于 侦破 命案 而 言 是 
极其 低 效 的 , 这 也 是 法 医 昆 虫 学 在 实际 办 案 中 很 难 
应 用 的 根本 原因 。 

随 着 DNA 检测 和 分 析 技 术 的 飞速 发 展 , 法 医 
昆虫 学 领域 逐渐 尝试 着 应 用 各 种 分 子 手段 对 法 医 昆 
虫 (尤其 是 嗜 尸 性 蝇 类 ) 进行 种 属 鉴定 ， 如 RFLP, 
AFLP 以 及 DNA 测序 等 (Sperling et al., 1994; 
Kakouli-Duarte et al., 2001; Amendt et al., 2004) 。 
近年 来 , 通过 测定 某 个 或 者 某 几 个 基因 的 部 分 序列 
来 进行 物种 鉴定 和 区 分 的 DNA. 条 形 码 ( DNA 
barcode) 技术 在 国外 得 到 了 迅速 的 发 展 和 应 用 
( Hebert et al., 2004b; 首 金 花 等 , 2004; Meyer and 
Paulay , 2005; Savolainen et al., 2005) 。 使 用 这 一 技 
术 可 以 让 非 专 业 人 员 在 很 短 时 间 内 准确 、 经 济 地 对 
目标 物种 进行 鉴定 和 区 分 , 大 大 提高 了 与 物种 鉴定 
有 关 的 学 科 如 分 子 生态 学 (Smith et al., 2005; 王 剑 
IE JU Zr f Be, 2007). P&P ^E 2] 5€ (Hebert et al., 


20042) 和 农业 科学 (Ward et al., 2005; Ez ids, 
2007) 等 的 应 用 效率 和 实用 价值 。 线 粒 体 DNA. 上 的 
细胞 色素 C 氧化 酶 亚 基 工 基因 (cytochrome oxidase 
subunit 1, COI) 的 100 ~650 bp 被 证 实 能 对 很 多 动 
物 物 种 进行 高 效 区 分 ( Hebert et al., 2003; Stoeckle, 
2003 ) 。 因 此 , 它 也 是 进行 法 医 昆 虫 学 DNA barcode 
研究 的 首选 分 子 材料 。 

本 研究 即 是 应 用 测序 技术 对 前 述 1993 年 所 确 
定 的 北京 地 区 篆 见 的 3 科 4 属 7 种 嗜 尸 性 蝇 类 线 粒 
体 基 因 COI 1 120 bp 的 DNA 序列 进行 检测 , 在 通 
过 序列 比较 对 有 关 个 体 进行 鉴定 和 区 分 的 同时 , 也 
为 后 续 DNA barcode 的 分 析 选 择 和 单位 点 分 型 研究 
提供 参考 信息 。 


1 材料 与 方法 


1.1 样本 采集 

哮 忆 性 蝇 类 样本 采集 目 2006 - 2007 年 北京 市 
案 发 现场 尸体 周围 具体 地 点 覆盖 顺义 (北京 东 
WB). EE CV). Bil Cg) o8 2g CAE) $8 HET] 
(中 ) 等 5 PEE AS RIZ I8] B DC S RE 8] EIS] NL HR 
样本 共计 77 3k, 样本 采集 后 浸泡 于 95% 的 乙醇 中 ， 
由 中 国 科学 院 劲 物 研究 所 张学忠 进行 形态 学 分 类 鉴 
E, 鉴定 结果 为 : AAWE k, KAEH 17 头 ， 
丝光 绿 蝇 19 3L, JEZER 7 3k, IX P NV RR HRS 1 
头 , 标 尾 别 麻 晶 6 头 和 家 蝇 26 k, 详情 参见 表 1。 
1.2 DNA 提取 

用 德国 QIAGEN 公司 的 QIAamp DNA Mini Kit 
AIR ER ERE N EIK mtDNA, BAIER 
方法 如 下 :(1) 样 本 前 期 处 理 : EE EXETEAS HAZ 
漂洗 4~5 次 ,吸水 纸 吸 干 后 ,分 别 用 消毒 剪刀 藤 去 
TEk HEKRA, 将 留 取 的 苍蝇 胸部 肌肉 放置 于 
HAF, 用 狠 子 夹 取 包 囊 于 锡 纸 中 , 迅速 用 研 钵 将 
锡 纸 内 肌肉 组 织 碾 成 粉末 , 分 别提 取 25 mg 样本 粉 
末 置 于 1.5 mL 离心 管 中 备 用 ; (2) 加 入 180 pL 
Buffer ATL; (3) 加 入 20 uL Proteinase K, 56€ 7K% 
消化 3 h, 每 小 时 振荡 2 ~3 次 ; (4) 轻 微 离 心 , 加 入 
200 uL AL, 充分 混 匀 , 70C 水 浴 10 min; (5) 加 入 
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200 uL (9696 ~100% ) ZB, 充分 混 匀 , 轻微 离心 ; 
(6) 转 移 混 合 物 到 QIAamp spin 中 , 8 000 rpm 离心 
1 min, 换 套 管 ; (7) 加 入 500 uL AW1, 8 000 rpm 
离心 1 min， 换 套 管 ; (8) 加 入 500 pL AWZ, 
14 000 rpm 离 心 3 min, 换 1.S mL 离心 管 ; (9) 加 入 
200 LÆ tak, 室温 静 置 min, 8 000 rpm 离心 
1 min, -20C 保存 备用 。 


1.3 引物 设计 

利用 Primer Premier 5. 0 软件 (Premier 公司 , 加 
拿 大 ) 设计 COI 引 物 , 扩 增 产物 覆盖 了 COI 基 因 
106 bp 到 1 326 bp 的 序列 。 具 体 引 物 序列 为 : 
Sense (F): 5'-CGAGCTGAATTAGGACATCCTGG-3' 23 bp 
Anti-sense (R): 5'-TGGGTAATCTGAATAACGTCGAGG-3' 24 bp 


R1 嗜 尸 性 蝇 类 样本 采集 信息 和 COI 基因 序列 的 GenBank 登录 号 


Table 1 Collecting data of sarcosaphagous flies in various localities in Beijing and 


GenBank accession numbers of COI sequence 


种 类 


Species 


丽 蝇 科 丝光 绿 蝇 
Calliphoridae 


Lucilia sericata 


红头 丽 量 Calliphora vicina 


大 头 金晶 
Chrysomyia megacephala 


PEREP PR FEJ RA 


Sarcophagidae Boettcherisca peregrina 


HESR MV. RA 


Parasarcophaga crassipalpis 


Parasarcophaga portschinskyi Rohdendorf 


BERI 2c 


Muscidae Musca domestica 


1.4 PCR 扩 增 和 测序 


PCR 反应 体系 为 25 uL, 具体 包含 : 2.5 uL 10 





























采集 地 点 采集 日 期 GenBank 登录 号 
Locality Date ( yy-mm-dd) GenBank accession no. 
顺义 Shunyi 06-07 -09 EU815010 - EU815013 
昌平 Changping 06-07 -11 EU815014 — EU815018 
房山 Fangshan 06-08-03 EU815019 — EU815023 
天 通 苑 Tiantongyuan 07-08-12 EU815024 
天 通 苑 Tiantongyuan 06-08-13 EU815025 - EU815026 
德 胜 门 Deshengmeng 07-08-22 EU815027 — EU815028 
天 通 苑 Tiantongyuan 07-09-27 EU814982 
顺义 Shunyi 06-07 -09 EU815035 — EU815036 
顺义 Shunyi 07-09-27 EU815037 
昌平 Changping 06-07 -11 EU815038 — EU815040 
房山 Fangshan 06-08-03 EU815041 — EU815044 
天 通 苑 Tiantongyuan 06-08-13 EU815045 — EU815047 
德 胜 门 Deshengmeng 07-08-22 EU815048 
顺义 Shunyi 06-07 -09 EU815029 
昌平 Changping 06-07 -11 EU815030 
房山 Fangshan 06-08-03 EU815031 
天 通 苑 Tiantongyuan 06-08-13 EU815032 — EU815034 
顺义 Shunyi 06-07-09 EU814975 — EU814976 
昌平 Changping 06-07-11 EUS814977 — EU814978 
房山 Fangshan 06-08-03 EU814979 — EU814980 
天 通 苑 Tiantongyuan 06-08-13 EU814981 
天 通 苑 Tiantongyuan 07-08-12 EU814983 
顺义 Shunyi 06-07 -09 EU814984 — EU814990 
昌平 Changping 06-07 -11 EU814991 — EU814995 
房山 Fangshan 06-08-03 EU814996 — EU814999 
天 通 苑 Tiantongyuan 07-08-12 EU815000 - EU815007 
德 胜 门 Deshengmeng 07-08-22 EU815008 — EU815009 


x buffer, 2. 5 uL dNTP (2.5 mmol/L), 1.6 uL 
Mg * (25 mmol/L) , 13.5 uL XXZ&;K, 3 uL DNA 模 


LT 
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板 , 0.3 uL Taq DNA 聚合 酶 (5 UL ,华美 公司 ) , 
正 反 引物 各 0.8 uL (5 nmol/L) 。 通 过 Perkin2Elmer 
9600 扩 增 仪 进行 PCR 扩 增 。 反 应 条 件 为 95 加热 
5 min, 然后 进行 35 次 PCR 循环 : 95% 变性 30 s, 
52*C 3B X 35 s, 72% 延伸 $$ s, 最 后 72% 延伸 7 min, 
PCR 产物 保存 于 4C,， 1% 33 BE P USE Rz FR. UK 23 A 
PCR 772], Goldview 核酸 染料 (EB 替代 品 ) 染色 。 
扩 增 产物 用 使 用 BigDye Terminator V3. 1 Cycle 
Sequencing 试剂 盒 进行 正 反 辐 测 序 反 应 ， 按 推荐 方 
法 操作 。 使 用 ABI 3130 (AB 公司 ,美国 ) 自动 测序 
仪 进行 测序 。 
1.5 序列 分 析 

利用 Lasergene v7. 2( Burland, 2000) 的 SeqMan f£ 
序 对 每 一 个 体 正 反 问 测序 的 结果 拼接 , 然后 使 用 
MegAlign 程序 比 对 并 手工 校对 同一 物种 的 多 个 样本 序 
列 , 最 后 通过 所 有 序列 的 比 对 和 比较 确定 每 一 个 体 可 
用 的 区 段 。 在 此 基础 之 上 , 使 用 75 bp 和 51 bp 的 滑动 
窗口 考察 所 有 序列 可 变 位 点 的 分 布 情况 , 窗口 和 和 窗口 
之 间 分 别 取 相应 窗口 长 度 的 L/3 进行 重合 。 使 用 
Mega4 软件 (Tamura et al., 2007) 采用 Kimura 2- 
parameter 模型 ( Kimura, 1980) 和 最 大 组 合 似 然 模 型 
( Maximum Composite Likelihood, MCL) ( Tamura et al., 
2004) 构建 序列 间 的 UPGMA (Sokal and Sneath, 1963) 
和 邻接 (Neighbor-Joining ) ( Saitou and Nei, 1987) 系统 
发 生 树 ， 同 时 采用 bootstrap $5 3s (Felsenstein, 1985) 
(1000 次 ) 计 算 系 统 发 生 树 每 一 分 文 的 支持 值 。 男 外 ， 
还 使 用 Mega4 计算 物种 内 和 物种 间 的 分 卜 情 况 , 计算 
采用 Kimura 2-parameter 模型 。 


2 fm 


2.1. 序列 结果 和 滑动 窗口 统计 

经 过 正 肥 向 拼接 、 同 物种 序列 比 对 和 手工 校正 
以 及 去 除 低 质 量 区 段 , 我 们 得 到 了 最 后 可 供 比 较 分 
析 的 高 质量 序列 77 条 , 序列 长 度 分 布 在 1 076 bp 
至 1 155 bp 之 间 , 平均 为 1 105 bp。 所 有 序列 提交 
至 GenBank 得 到 相应 的 编写 ( 表 1)。 后 面 的 分 析 以 
多 序列 比 对 的 1 076 bp 公共 区 段 为 参照 。 

通过 滑动 窗口 考察 序列 可 变 位 点 分 布 的 结果 见 
图 1。 从 图 中 我 们 可 以 看 到 , 序列 间 的 差异 位 点 在 
整个 研究 区 段 内 均 有 分 布 且 密 度 差别 不 大 , 位 扣 数 
随 窗 口中 心 的 移动 有 轻微 下 降 趋 势 。500 bp 附近 
可 变 位 点 数 最 少 , 而 600 ~ 650 bp 附近 可 变 位 点 数 
最 多 。 所 以 , 相 比 使 用 部 分 序列 进行 比较 , 使 用 序 


列 整个 区 段 能 够 包含 更 丰富 的 信息 。 


一 一 75 bp 滑动 窗口 75 bp sliding window 


3Ó[ — ---- 51 bp 滑动 窗口 51 bp sliding window 
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可 变 位 点 个 数 
Number of variable sites 





0 200 | 400 600 


800 1000 1200 
窗口 中 心 位 置 bp) 


Center of windows 
图 1 序列 可 变 位 点 的 滑动 窗口 分 析 
Fig. 1 Sliding window analysis of variable sites 
among all sequences 
相 邻 窗口 互相 覆盖 1/3 窗口 长 度 。Adjacent windows have an overlap 
of 1/3 window span. 


2.2. 系统 发 生 树 构建 

图 2 是 77 条 序列 间 的 UPGMA 系统 发 生 树 。 
图 中 分 属 各 物种 的 个 体 的 聚集 趋势 十 分 显著 , 相应 
分 支 的 支持 值 均 超 过 99% , 表明 同一 物种 的 个 体 可 
以 同 其 他 物种 的 个 体 很 好 地 区 分 开 。 在 种 间 水 平 ， 
家 晶 同 丽 蝇 科 的 3 个 物种 以 及 腑 蝇 科 的 3 个 物种 间 
的 分 文 文 持 值 分 别 为 95% 和 92% , Tii WE EPA R 
量 科 内 部 的 分 支 支 持 值 只 有 55% 和 76% 。 其 中 肥 
须 亚 肪 晶 和 急 钧 亚麻 晶 同 属 亚麻 晶 属 , 但 别 麻 晶 属 
的 标 尾 别 麻 晶 和 急 钩 亚麻 晶 表 现 出 更 近 的 距离 ,， 且 
相应 分 文 仅 有 76% 的 支持 值 。 

我 们 将 样本 序列 与 一 段 长 278 bp 的 COI 序列 
AY818110 (ZAKIS, 2005) 的 重合 区 有 段 做 了 类 似 
的 系统 发 生 分 析 。 该 重合 区 段 的 长 度 只 有 186 bp, 
但 是 得 到 的 系统 发 生 树 在 物种 水 平 的 分 文 模式 与 
1 076 bp 全 序列 的 结果 一 致 , 且 各 物种 内 的 序列 仍 
以 99% 的 分 文 文 持 值 聚集 在 一 起 (数据 未 展示 ) o 

另外 我 们 选用 序列 的 最 大 组 合 似 然 模 型 对 
1 076 bp 和 186 bp 序列 构建 UPCMA 发 生 树 , 结果 得 
到 同样 的 分 文 模 式 ， 同 一 物种 的 序列 以 100% 的 支持 
值 聚 集成 族 ， 而 种 间 分 支 的 支持 值 存在 略微 的 变化 
(图 2 以 及 未 展示 数据 ) 。 当 采用 邻接 方法 构建 发 生 
树 时 , 大 头 金晶 对 应 的 复 和 家 蝇 篮 接合 在 一 起 (数据 
未 展示 ) 。 这 一 结果 和 分 类 学 结论 存在 分 歧 , 不 过 上 
述 接 合 分 文 的 支持 值 为 90% (对 应 1 076 bp 结果 )， 
IRF UPGMA 发 生 树 中 大 头 金晶 和 丝光 绿 蝇 艇 相 接 
44] x HJ 9596 支持 值 。 改 用 Kimura 2-parameter 模型 
同样 也 不 影响 NJ 建树 方法 的 结果 。 
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图 2 1076 bp 序列 给 出 的 个 体 间 UPGMA 系统 发 生 树 
Fig. 2  Unrooted UPGMA tree constructed under Kimura 2-parameter model with 1 076 bp sequences 
Jia: 家 晶 Musca domestica; D ( DT): AK 3k 4 W& Chrysomyia megacephala; FY: HE Wi WV Jk WB Parasarcophaga crassipalpis; JG: 4 £5 WV. Wk Wü 
Parasarcophaga portschinskyi Rohdendorf; ZB: fx E ppkt Boettcherisca peregrine; H: ZL3k Wii Calliphora vicina; S: 2£J65&W& Lucilia sericata. 图 中 
的 数值 为 1 000 次 bootstrap 检验 的 支持 值 , 同 物种 个 体 序列 聚集 的 支持 值 均 达到 99% 。 种 内 序列 分 支 的 支持 值 未 予 显示 。Sequences of same 
species are highly clustered with a support value of 99%. The number for each interior branch is the bootstrap value, and values within species are not 


shown. 
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2.3 序列 分 政 统 计 

物种 内 部 和 物种 间 的 序列 分 监 统计 见 表 2。 从 
表 中 可 以 看 到 , 序列 的 种 内 差异 为 0.0016(S) 到 
0.0069(ZB) 不 等 , 均 未 超过 1%。 种 间 差 异 净值 最 
小 出 现在 ZB 和 JG 之 间 , 为 0.0774, 最 大 值 出 现在 
JIA 和 FY 之 间 , 达到 了 0. 1485。 此 结果 和 系统 发 
生 分 析 的 结论 相 一 致 , 种 内 序列 差异 较 小 因而 会 在 
所 有 序列 的 分 析 中 紧密 聚集 在 一 起 。 

TFEZRGK SE RU QI FE 13-8] £T 3E WI BRURI 3.44] MI JPRÉER EIS] 
DNA FI 1 条 。 为 了 在 种 内 序列 分 时 统计 上 不 
影响 和 其 他 物种 的 比较 , 我 们 将 这 两 个 DNA. 序列 
在 GenBank 做 了 blastn 比 对 。 其 中 和 红头 丽 蝇 序列 
H1 相似 性 最 低 的 同 物种 序列 为 DQ345096 ， 达 到 
99. 596 ; 序列 相似 性 最 高 的 不 同 物种 序列 为 
DQ345095 ( S iji Calliphora pattoni) , 仅 为 96.49% 。 
PEIRE EIS] 37 20] JG6 相似 度 最 高 的 不 同 物种 序 
列 是 EF405937( 桨 亚麻 蝇 Sarcophaga dux) ， 相 似 度 
为 94.4% ; 同 物种 序列 未 找到 。 这 些 结果 表明 , R 
AZL3E WE ERCRT Za £9 MT PREIS] FE 91] 43-90] FUE 1 条 , 但 
是 我 们 得 到 的 序列 信息 足以 将 它们 同 其 他 物种 区 
分 开 。 

R2 种 内 和 种 间 序 列 分 歧 情 况 
Table 2 Sequence divergence between and within species 
FY H JG Jia S ZB D 

FY 0.0027 0.1232 0.0867 0.1485 0.1113 0.099 0.1131 

H 0.1245 NA 0.1304 0.1341 0.0846 0.1326 0.1028 
JG | 0.0881 0.1304 NA 0.1408 0.1033 0.0774 0.1211 

Jia 0.1509 0.1351 0.1418 0.002 0.119 0.1428 0.1069 

S | 0.1134 0.0854 0.1041 0.1208 0.0016 0.1158 0.0848 

ZB 0.1037 0.136 0.0808 0.1472 0.1201 0.0069 0.1291 

D 0.1157 0.1041 0.1224 0.1092 0.0869 0.1339 0.0026 
种 名 代码 同 图 2。 对 角 线 为 种 内 序列 差异 , 下 三 角 为 种 间 差 异 , 上 
三 角 为 除开 种 内 差异 的 种 间 差 异 净值 。NA 表示 仅 一 条 序列 。 
Species code is the same as Fig. 2. Diagonal, lower triangle and upper 
triangle values are for within-species, between-species and net between- 
species divergence, respectively. NA means just one sequence and no 


comparison. 
3 讨论 


通过 进行 线粒体 DNA E COI 基因 1 076 bp JE 
列 的 测定 和 比较 分 析 , 不 同 地 点 采集 的 不 同 种 类 嗜 
尸 性 晶 类 个 体 得 到 了 成 功 的 区 分 , 序列 信息 的 聚集 
结果 和 形态 学 鉴定 的 种 类 归属 完全 一 致 ， 因此 本 研 
究 得 到 的 各 物种 的 序列 完全 可 以 作为 DNA. barcode 


应 用 于 实际 案件 的 调查 。 后 续 实 际 应 用 中 , 对 案件 
现场 采集 的 嗜 尸 性 蝇 类 蛆虫 或 者 成 虫 提取 线粒体 
DNA 后 进行 测序 , 并 与 我 们 的 barcode 序列 进行 比 
对 ， 即 可 根据 序列 的 分 支 归 属 确定 其 物种 。 这 种 方 
法 的 推广 在 保证 物种 鉴定 准确 性 的 同时 , 将 能 大 大 
降低 嗜 尸 性 蝇 类 在 法 医 鉴定 中 的 应 用 难度 。 

另 一 方面 , 滑动 窗口 分 析 显 示 可 变 位 点 信息 基 
本 遍布 CO 基因 的 整个 区 段 , 同时 186 bp 序列 比 
对 分 析 与 全 序列 给 出 的 区 分 一 致 (化 继 峰 等 ， 
2005), 4 F COIL 基因 的 可 变 位 点 分 布 特征 ， 
650 bp 甚至 更 短 的 片段 序列 即 可 达到 类 似 使 用 全 序 
列 时 的 区 分 效果 (Hajibabaei et al., 2006b) ， 所 以 不 
同 引 物 设计 得 到 的 部 分 区 段 不 论 长 短 , 均 可 以 和 本 
研究 给 出 的 长 序列 进行 比较 ( 即 读 取 barcode 过 
T£), ， 从 而 实现 有 关 物 种 的 区 分 和 鉴定 。 

在 使 用 全 序列 或 部 分 序列 构建 的 系统 发 生 树 
中 , 同一 物种 个 体 序列 聚集 的 同时 也 有 部 分 物种 整 
体 出 现在 和 分 类 学 结果 不 一 致 的 位 置 。 例 如 ,同属 
XV PR RES EI] i 4] MI. JPR BER. 3 HS 20 NIV. JB B < [R] A ER S Jc 
m EA EJE B0 3.4 MI. JR B 3 ps EE PE [|] EA] EB S 
表现 的 更 远 。 所 以 虽然 应 用 了 很 多 系统 发 生 的 工 
H, 但 在 涉及 进化 树 的 这 度 构建 时 ，parcode 并 不 能 
代 符 传统 的 系统 发 生 分 析 ,， 分 析 barcode 序列 构建 
的 系统 发 生 树 不 能 下 接 视 为 相应 物种 在 进化 上 的 系 
统 发 生 关 系 , 这 和 之 前 的 很 多 研究 是 一 致 的 (Russo 
et al., 1996; Zardoya and Meyer, 1996; Hajibabaei et 
al., 2006a) , 

各 种 DNA 分 子 标记 , 包括 STR Y 染色 体 
SNP 等 等 , 被 广泛 应 用 于 亲子 鉴定 和 个 体 识别 等 领 
域 ( 李 备 华 等 , 2001; 李 生 斌 等 , 2001)。 从 本 质 上 
KH, 这 都 是 利用 了 亲缘 关系 远近 决定 遗传 信息 差 
腊 的 水 平 。 本 研究 所 获得 的 约 1 100 bp 的 序列 也 可 
以 为 筛选 用 于 嗜 尸 性 晶 类 物种 区 分 的 线粒体 DNA 
差异 位 点 提供 帮助 。 通 过 考察 所 有 哮 尸 性 蝇 类 每 一 
物种 多 个 独立 个 体 的 序列 , 我 们 很 容易 确定 在 线 粒 
体 DNA 某 一 区 段 某 一 物种 所 特有 的 突变 位 点 的 组 
合 ,， 从 而 根据 对 这 些 组 合 位 点 的 分 型 达到 实现 区 分 
相应 蝇 类 物种 的 目的 。 这 样 的 研究 在 其 他 如 细菌 、 
水 称 等 领域 已 经 有 一 些 研 究 工 作 (Troesch et al., 
1999 ; GAMA, 2003; 2E XC ARS, 2008), ,而 在 
法 医 昆 虫 学 的 研究 中 此 方法 具有 同样 重要 的 开发 价 
值 。 对 嗜 尸 性 蝇 类 种 类 和 时 间 分 布 的 进一步 研究 、 
新 的 基因 分 型 技术 的 发 展 以 及 有 关 物 种 线粒体 
DNA COIL 基因 序列 的 进一步 积累 ,最 终 将 帮助 我 们 
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选择 最 佳 的 突变 位 点 组 合并 由 此 实现 更 加 便利 的 嗜 
己 性 蝇 类 物种 的 鉴定 。 
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